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Abstract 


Prior information for the stationary sea surface topography (SST) may be 
needed in altimetric solutions that intend to simultaneously improve the 
gravity field and determine the SST. For this purpose the oceanographically 
derived SST estimates are represented by a spherical harmonic expansion. 
The spherical harmonic coefficients are computed from a least squares 
adjustment of the data covering the majority of the oceanic regions of the 
world. Several tests are made to determine the optimum maximum degree of 
solution and the best configuration of the geometry of the data in order to 
obtain a solution that fits the data and also provides a good spectral 
representation of the SST. 
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Introduction. 


One of the main problems in marine geodesy is the determination of the 
geoid from the sea surface heights that are computed from satellite altimetric 
observations. In order to do that independent estimates of the stationary and 
time variable sea surface topography (SST) need to be provided. Time 
variations of SST can be computed in a rather straightforward way by 
analysing the overlapping tracks during the repeat era of a satellite altimetric 
mission (Cheney et al., 1983). These variations can then be removed from the 
sea surface heights. Alternatively they can be considered as noise during the 
processing of altimeter data with traditional crossover techniques and 
therefore be filtered out from the sea surface heights (e.g. Rowlands, 1981). 
Estimates of the stationary SST, on the other hand, that has an expected total 
variation of about 2 meters, can be provided by oceanographic methods (e.g. 
Levitus, 1982) using observations of ocean temperature, salinity and oxygen 
content and imposing geostrophic conditions in solving the equations of motion 
for the oceans. 

Stationary SST can also be computed by traditional geodetic techniques. 
In one such determination (Engelis, 1985) geoid undulations computed from a 
low degree satellite derived gravity field, are subtracted from the sea surface 
heights computed from altimetry. Because of the errors in the determination 
of the satellite derived gravity field only the long wavelength part of these 
reduced sea surface heights can represent the stationary SST with some 
degree of confidence. In order to determine this long wavelength SST an 
harmonic analysis of the reduced heights is required. Then a low pass 
filtering is performed to retain only the low degree coefficients that have a 
favourable signal to noise ratio. 

Recently there have been attempts to simultaneously determine the long 
wavelength SST and improve the gravity field of the earth. One such method 
incorporates a low degree (n max = 10) spherical harmonic model into a general 
dynamic solution for a low degree (n max = 50) gravity field of the earth, 
currently being attempted at NASA Goddard Space Flight Center. In such a 
solution observations to geodetic satellites, altimeter data and terrestrial 
gravity anomalies are used. In an alternative method being proposed by 
Engelis (1987), altimeter observations are used in a combined solution to 
reduce the radial orbit error, improve the geoid and determine the stationary 
SST. 

The effectiveness of all geodetic techniques to determine the SST is 
subject to the accurate representation of the spectral content of the SST 
itself, or, in other words, subject to the correct implementation of the 
spherical harmonic models that are used. Particularly for the last method, 
where apriori SST information by wavelength (i.e. degree variances) is 
needed, it is important that some estimates of the spectral behavior of SST is 
available. Such estimates have been provided in the past by Engelis (1985) 
who has used the Levitus data to estimate harmonic coefficients and their 
degree variances. In that estimation the orthogonality principle assuming data 
all over the sphere was used. In the present investigation this determination 
is reexamined since the SST is not a complete function on the sphere but is 
only defined in oceanic areas. Therefore, consideration of the SST as a global 
function introduces problems, the most important of which is that, the 
resulting harmonic coefficients have a lower power since they are forced to 
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also fit the land regions which are filled with SST values that are traditionally 
considered to be zeroes. In the present analysis the Levitus data, only in the 
ocean areas, will be used to determine spherical harmonic coefficients. These 
coefficients can provide, hopefully, the best reference values needed for the 
combined solutions. A more general purpose of this analysis is to examine the 
problems involved with the spectral decomposition of data that are not globally 
distributed on the sphere. 


The Levitus SST. 


In oceanographic methods to determine the stationary SST, temperature, 
salinity and disolved oxygen content of the oceans are used to determine 
pressure and water density, which in turn are used to solve the geostrophic 
equations of motion in the oceans. In order to do that a reference 
equipotential surface (surface of no motion), which ideally can be the geoid, 
must be used. Due to the inability to use an estimate of the geoid, a deep 
surface is defined to be a surface of no motion. Then solution of the 
equations of motion provides the mean annual anomaly of the geopotential 
thickness of the layer between that surface and the ocean surface. 

The most recent determination of such a dynamic SST is made by Levitus 
(1982) who used data from the National Oceanographic Data Center. A first 
analysis of the data, made by Levitus, indicated several problems. The most 
important problem was regional biases in the data and lack of data in 
extended regions. Moreover there were temporal representation problems in 
the data since observations were not synoptic but scattered with respect to 
time (with the exception of a few limited areas) and so, the results cannot in a 
strict sense be considered a true long term average. 

After the initial stability and statistical checks to eliminate spurious 
observations, averages of data in l’xl* blocks were created. In order to 
overcome biases and lack of data a smoothing operation was performed. This 
smoothing consisted of a weighted average operation in which a Gaussian type 
filter was used as a weighting operator. The radii of the Gaussian filter 
ranged from 1540 km to 770 km depending on the region. As a result, any 
signal with wavelengths less than 800 km was eliminated while wavelengths 
between 800 km and 3000 km were affected with changes in the amplitudes of 
the signal. For example wavelengths of 1000 km had a reduction in amplitude 
of at least 50% (Levitus, 1982, Figure 11). The minimum wavelengths of 800 km 
roughly correspond to a maximum degree 20 in a spherical harmonic expansion. 
Levitus considers that the resulting large scale features are representative of 
the real ocean, although it is expected that some local differences can occur 
because of interannual variability. This smoothed data set was used to 
compute the annual mean anomaly of geopotential thickness of several layers 
corresponding to different deep surfaces considered to be surfaces of no 
motion. 

The dynamic topography that is used in the present analysis is the one 
with respect to a 2250 db surface. This data set consists of 33856 l*xl* mean 
values in the ocean areas of the world. Their spatial distribution is shown in 
Figure 1 of Appendix A. A first analysis of this data set has indicated some 
outliers in the west equatorial region of the Pacific ocean. After these values 
are rejected, the weighted mean value of the SST set is computed to be 2.02 
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meters. This mean is removed, since for geodetic applications any such terms 
are absorbed by the mean earth ellipsoid that is used to reference the geoid. 
Analysis of the centered quantities indicates that the SST estimates range 
from -1.40 meters at the southernmost latitudes close to Antartica, to about 1 
meter in the northern Pacific. A notable exception to the above range has 
been found in the Mediterranean sea where there is a sparse coverage of 82 
l*xl* values all of them having magnitudes smaller than -3 meters. 
Additionally in the oceanic regions above 70* as well as in the region between 
Greenland and Scandinavia there are 3134 SST estimates with a mean value of 
-1.80 meters and a very small, almost latitudinal variation. These values differ 
from neighboring values in the North Atlantic region as well as values in the 
Northern Pacific by as much as 1 meter. 

This original data set with the mean value removed is the one that was 
used in Engelis (1985). In expanding the SST into spherical harmonics one is 
interested in determining a set of coefficients that best represent the oceanic 
regions within latitude limits that are also attained by satellite altimetry. 
Therefore, it was decided (Rapp, 1985, Appendix B) to reject all the 3134 
values in the northernmost latitudes. Furthermore, the Levitus estimates in 
the Mediterrannean sea were replaced by estimates from a map by Lisitzin 
(Lisitzin 1974, p.153) that was given with respect to a 4000 db surface. In 
order to put the Lisitzin estimates into the same reference system as the 
Levitus set, their mean value was removed and the global mean value of the 
Levitus set was added. 

In the present investigation four sets of SST estimates are used in order 
to see what is the effect of particular regions on the global analysis of SST. 
These sets are the following. 

1. Original Levitus data set centered around its mean value of 2.02 meters 
and containing 33856 estimates with an RMS value of 80.3 cm (SET1). 

2. Same as SET1 but augmented by the Lisitzin estimates in the 
Mediterrannean sea, centered around its mean value of 2.02 meters and 
containing 34056 estimates with an RMS value of 78.8 cm (SET2). 

3. Same as SET2 but without the 3134 estimates in the northernmost 
latitudes. There is a total of 30922 estimates with their mean value of 
2.01 meters removed and an RMS value of 62.4 cm (SET3). 

4. Same as SET3 but with no data in the Mediterrannean sea. There are 
30640 estimates with their mean value of 2.03 meters removed and an 
RMS value of 62.4 cm (SET4). 

The spatial distribution of all four data sets is shown in Figures 1-4 of 
Appendix A. 
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Harmonic Analysis of a Function on a Sphere. 


A square integrable analytical function f(4»,X) defined on a unit sphere 
- 7 T/ 2 < 4 ><rr /2 and can be expanded in a series of surface spherical 

harmonics 


f(*,X) 



where 


C im 

if 

a-0 

Sim 

if 

a ; =l 


and 

P^ m (sin<t>)cosmX if a=0 

P^ m (sin4>)sinmX if a=l 



( 1 ) 


( 2 ) 


(3) 


are the fully normalized spherical harmonic coefficients of the 
function f (4>,X) and Pi m (sin4>) are the fully normalized associated Legendre 
functions of the first kind, such that 


1_ f =a rb , (1 if 1=1, m=j, a=b 

4tt J 'j Q lo otherwise 


(4) 


Using the orthogonality principle that is expressed by equation (4) the 
spherical harmonic coefficients can be derived from 

Ci m = ^ | *<♦,*) Yl rn (4>,X)d < x (5) 


These coefficients are independent and they form a basis of the function on 
the sphere. In terms of linear algebra, equation (1) is the spectral 
decomposition of the function f(<t>,X), are the eigenfunctions and Cj> m are 

the eigenvalues of the function. 

When discrete point realizations of the function are given on the sphere 
in an equiangular gridded form, with grid intervals A4>, AX along latitude and 
longitude respectively, then the upper limit of equation (1) changes from 
infinity to some maximum degree i max that corresponds to the Nyquist 
frequency of the sampled set and is equal to tt/AX. Furthermore, equation (5) 
becomes 


L-i u-i 


°5.= t xj>«u 


where <Tjj is the finite surface element and is equal to 


( 6 ) 


= AX(sin(4>j+A<l>)-sin<l> * )' 


(7) 
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In equation (6), L is the number of gridded samples along a parallel and is 
equal to ^ max . The orthogonality principle is still valid and independent 
coefficients can be obtained up to ^ max . These coefficients are insensitive to 
the maximum degree of solution as long as this is smaller than i max since 
folding of frequencies is prevented by the orthogonality of the harmonics. 
Attempts to determine coefficients of higher degrees result into aliased 
estimates that are fully correlated with coefficients of degrees lower than 

i 

max* 

When the function f (4»,X) is sampled in the form of equiangular area means, 
things become a little more complicated. Equation (1) becomes 


1 Imax £ i — a 

f(*i,A<) = — E E E c!_ 

1 J J=o m=o a =o 


f Y| m d<r 


« J 


or, in an expanded form 


( 8 ) 





£max 

ii. 


£ 

E 


m-o 


(C £m IC m +S £m IS m )IP £m 


(9) 


where IC m , IS m , IPi m are the integrated cosines and sines of order m and the 
integrated associated Legendre functions, respectively. The maximum degree 
of expansion is again specified by the relationship i max = 7r/AA where AA is 
now the size of the block. 


The problem in expanding area means into spherical harmonics is that the 
orthogonality principle is not valid anymore. As a matter of fact deviations 
from the orthogonality principle are a function of the block size. These 
deviations tend to zero as the block size tends to zero. For a given block 
size deviations from the orthogonality principle increase with increasing 
degree. Approximations to the orthogonality principle and therefore estimates 
of the spherical harmonic coefficients can be obtained by introducing the so 
called desmoothing factors (Colombo, 1981, Rapp, 1986). Then equation (5) 
becomes 


c u = 4^- V 2 V f(-Mj) J f \j)d* ao) 

where qi can be either the Pellinen operator or an optimum quadrature weight 
defined by Colombo (1981). With the incorporation of the desmoothing factor 
errors due to the finite size of the block and errors arising from the 
dampening of higher frequencies because of the averaging to generate the 
mean values, are reduced. 

Spherical harmonic coefficients computed using equation (10) are not 
independent any more due to the approximations involved in the equation 
itself. Their correlations though are minimal. For the efficient evaluation of 
(10) optimum procedures have been developed and can be found in Rapp 
(1986). 
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Spherical harmonic coefficients representing the square integrable 
analytical function f(4>,A) can also be computed by a least squares adjustment 
of the sampled values of the function. Equations (1) or (8) can be used for 
that purpose as observation equations depending on whether the sampling is 
in a form of point values or area means. In a least squares solution 
individual data errors can be taken into account for weithting purposes, or all 
weights may be set equal to one. In the latter case and for point values, the 
formulation and the results are identical to the ones of equation (6) since the 
orthogonality principle is valid and so the normal matrix becomes diagonal. 
When mean values are used, with equal weights, the normal matrix is 
dominantly diagonal and the results and formulation are very similiar to the 
ones of equation (10) but not identical. In either formulation the coefficients 
can be considered as independent and insensitive to the maximum degree of 
solution. Furthermore, when the data sampling implies a Nyquist frequency 
that is greater that the maximum frequency existing in the function data 
(when the function is band limited), all coefficients corresponding to 
frequencies between these two frequencies will be effectively zero. Obviously, 
after the coefficients are estimated, a subset of them can be used to compute 
a long wavelength approximation of the function (i.e. a low pass filtering) or, 
equally well, high degree variations (i.e. a high pass filtering). Finally the 
spectrum of the function can be computed by computing the degree variances 
as follows 


2 



£ 

i 


m-o 



+ s’ ) 

x nrr 


(ii) 


When the function f(<f>,A) is not globally defined on the sphere but only on 
a portion of it, then it is not analytical on the sphere anymore since its 
spatial derivatives are discontinuous on the boundaries. Consequently it 
cannot be expanded into a series of surface spherical harmonics that are 
orthogonal, since the orthogonality principle is not valid anymore. Ideally, 
such a function can be spectrally decomposed, if one is able to compute the 
eigenfunctions and eigenvalues of the function over the domain where the 
function is defined. When the boundaries though are variable, as it happens 
in most geodetic applications for which data are not sampled globally (e.g. 
satellite altimetry ) such a computation becomes nearly impossible. 

Furthermore, it has been traditional for geodesists to work with spherical 
harmonics since all the geodetic quantities of interest can be analysed into 
spherical harmonics. 

Based on the above, one can expand f($,A) into spherical harmonics, only 
with the understanding that such an expansion is nothing else than a 

polynomial fit, and that the only property of the harmonic coefficients is that 
they reproduce the function. To compute the coefficients of equations (1) or 
(9), depending on the nature of the available sampling, one cannot use 
equations (6) or (10) anymore. The only way these coefficients can be 

computed is through a least squares fit of the data using (1) or (9) as 
observation equations. A similiar situation arises when a function that is 
globally defined on the sphere is not globally sampled (e.g. terrestrial gravity 
anomalies). 

There are several issues that need to be addressed during such a 

computation. First of all and most important, the coefficients themselves are 
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not independent but correlated. This is expected since the normal matrix is 
not diagonal or even dominantly diagonal. Therefore, no long wavelength 
approximation to the function can be computed and equation (11) cannot really 
be used to compute the spectrum. Furthermore, the Nyquist frequency cannot 
be defined, based on the sampling of the function, in the sense i max = tt/AA 
since the discontinuities at the boundaries introduce artificial energies at 
frequencies that are functions of the variability of the boundaries. This 
effect is the so called Gibbs phenomenon, which is very well known in the 
theory of Fourier series. So even if the function is band limited with a 
maximum frequency i less than the maximum frequency i max that is implied by 
the data, if one tries to solve up to i max then the coefficients between £ and 
i max not b e zero but will have substantial magnitudes. Furthermore, 

considerable folding of these frequencies occurs so that to also affect the 
coefficients of degrees lower than i. As it will be seen later on, the greater 
the maximum degree of the solution is, the larger the folding and therefore 
the larger the coefficients become. Of course the correlations between these 
coefficients aldo increase in such a way that, when they are used to 
reproduce the function, they provide a perfect fit. 

From this discussion it is obvious that the computed coefficients are not 
unique but they depend on the maximum degree of solution. So, if a band 
limited function is analyzed, or smoothing has been applied to the data, it is 
important that some prior information about the effective minimum wavelengths 
(and therefore the maximum degree of solution) be available. 

In analyzing an incomplete data set on the sphere, either because of 
definition or because of incomplete sampling, the geometry of the data 
distribution is very important, and is basically the factor for the deviation of 
the harmonics from orthogonality. As will be seen later on, a very small 
change in the data distribution can result in substantial changes in the 
computed coefficients even when the same set of coefficients is solved for (i.e. 
same maximum degree of solution). Unfortunately there is no known method 
that can estimate what the effect of a gap in the data on the computed 
coefficients can be. The only information that can be obtained is how a 
particular gap of data affects the correlations of the coefficients. This can be 
done by simply computing the normal matrix for the gap. Doing that 
repeatedly, for several gaps, one can possibly identify areas to which a 
solution is very sensitive and areas that contribute a little to the solution. 


Determination of the SST Spherical Harmonic Coefficients. 

In determining the Levitus SST harmonic coefficients Engelis ( 1985 ) has 
adopted the definition that the SST is a global function that takes values both 
on land and oceans, with the values on land being identically equal to zero. 
In that definition it was required that the mean value of the SST is zero as 
sampled globally in the oceanic areas of the world. The SST estimates that 
were used in such an ivestigation were the original Levitus values that 
correspond to SET1. These estimates were considered to be point values. 
Then, spherical harmonic coefficients were computed using equation (6). The 
same analysis was repeated in Rapp (1985) where though SET3 was used and 
the values were considered to be mean values. Equation ( 10 ) was then used 
to compute the SST coefficients. In both solutions the results were practically 
insensitive to the maximum degree of expansion. The difference between the 
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two sets of coefficients was found to be on the order of millimeters for most 
of the coefficients with the exception of the first and second degree terms 
that differ by a couple of centimeters. Furthermore the cumulative power 
computed by either solution (up to 4 max = 36) was on the order of 40 cm 
which is substantially lower than the RMS values of the Levitus sets SET1 and 
SET3 that are on the order of 80.3 and 62.4 cm respectively. This substantial 
reduction was due to the smoothing of the derived fields that were forced to 
also fit the zeroes on land. For the same reason the RMS fit of these fields 
to the data was only on the order of 10-20 cm. 

The same type of solution is presently repeated by applying a least 
squares fit to the Levitus data set augmented by zeroes on land. Several 
solutions up to different maximum harmonic degrees (6,10,20,36) have been 
made for all the four sets of SST described previously. In all the solutions 
equal weights were used. Equation (9) was used to form the observation 
equations. The conclusions that can be drawn from the adjustment results are 
identical to the ones drawn from the quadratures solutions described above. 
Again there is ho variation of the coefficients among solutions with different 
£ max f° r any of the data sets. Furthermore solutions using SET1 and SET2 as 
well as solutions using SET3 and SET4 are almost identical (differences on the 
order of millimeters). This level of differences was also found in comparing 
these solutions with the corresponding quadratures solutions. In any of these 
least squares fits the standard deviation of each individual coefficient is 2 cm 
while the correlations were found to be indeed negligible. Any solution up to 
degree 36 has an RMS power of 40 cm and an RMS fit to the Levitus data of 
11 cm. Solutions up to degree 10, on the other hand, have an RMS fit of 20 
cm and an RMS power of 37 cm. 

Solutions for the SST harmonic coefficients have also been made by 
considering the SST to be defined only in the oceanic regions. In order to 
examine the sensitivity of the coefficients to the choice of the maximum 
harmonic degree and the geometry of the data distribution, solutions up to 8, 
10, 12, 20 and 36 were made using all the four data sets. In establishing 
which set of coefficients is the best, a combination of the following criteria 
was used. 

1. RMS fit to the Levitus data to be small. 

2. RMS value implied by the cumulative degree variances to be close to 
the RMS value of the Levitus SST. 

3. Zero degree coefficient to be very close to zero. 

4. Estimated standard deviations of the coefficients to be small. 

5. The condition number of the normal matrix to be small. 

This set of criteria was devised as an empirical alternative way to choose the 
set with the smallest possible correlations since examining the normal matrices 
for all these sets was found to be impractical. In particular, the combination 
of the first and second criteria can give a first indication about the level of 
correlations of coefficients. The third criterion can support the first two. 
Indeed, by definition the zero degree term has to be zero, so any deviation 
from zero indicates the level of correlation of the zero degree term with the 
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other coefficients. Finally the fourth criterion is useless by itself but if 
combined with the other criteria and particularly with the last one it can give 
one further input on the choice of the best set of coefficients. 

Examining all the different solutions it was found that the solution up to 
harmonic degree 10 using SET3 performs the best, since it better satisfies all 
the criteria. More specifically the RMS discrepancy from SET3 is 7 cm, the 
RMS value implied by the degree variances is 62.8 cm (closest to the RMS 
value of SST which is 62.4 cm) and the zero degree term is -4 cm. The 
standard deviations of the coefficients (not scaled by the aposteriori variance 
of unit weight) range from 2 to 30 cm and the condition number is on the 
order of 5000. These accuracy estimates are poorer than the corresponding 
ones from solutions to a smaller maximum degree (i.e. 8) but only marginally. 
On the contrary they are much smaller than the standard deviations of higher 
degree solutions. 

The performance of the solutions up to harmonic degree 10 using SET1 
and SET2 is almost comparable. Indeed the first two criteria are satisfied 
since the RMS fits to the corresponding data are on the order of 9 cm and 
their cumulative power reaches 80 cm. The other criteria though are not 
satisfied as well since the zero degree terms of the two solutions are -25 cm 
and -20 cm respectively and the standard deviations of all the coefficients are 
somehow larger than those of the SET3 solution. The reason for this is that 
the systematically lower values at latitudes greater than 70* introduce a step 
discontinuity in the data. This discontinuity affects primarily the zonals and 
particularly the second degree zonal that is on the order of -55 cm as 
compared to -28 cm of the SET3 solution. The coefficients computed from 
SET4 perform the worse. Indeed some of these coefficients have magnitudes 
greater than 50 cm giving a cumulative power in excess of 1 meter. Still, the 
RMS discrepancy from the SET4 data set is 7 cm. This indicates that there 
are high correlations between the coefficients and that the solution is 
extremely sensitive to the geometry of the data distribution. Indeed a small 
change of SET3 (removal of 282 values in the Mediterranean) makes the 
solution (i.e. SET4) completely unreliable. 

For all the data sets, any solution with a maximum harmonic degree 
greater than 10 gives coefficients that are completely unreliable, since they 
reach magnitudes close to 1 meter and standard deviations of several meters, 
while the condition numbers increase dramatically. In all these solutions the 
RMS discrepancies from the Levitus data reduce with increasing i max to reach 
a value of 8 mm for a solution up to degree 36. For this particular solution 
the cumulative power is on the order of 8 meters, the standard deviations are 
on the order of 10-15 meters and the condition number is on the order of 10 ®. 
Even a solution up to degree 12, although it gives an RMS fit of 6 cm to the 
Levitus data, it has a cumulative power of 1.35 m and individual standard 
deviations on the order of 50 cm. 

The coefficients of the ocean solution up to degree 10 and based on SET3 
are shown in Table 1 together with their by degree and cumulative amplitudes. 
Additional sets that are given for comparison are the ones of the global 
solution up to degree 10 (Table 2), and the ocean solution up to degree 36 
(Table 3), always based on SET3. All these Tables can be found in Appendix 
B. ' Comparing the degree amplitudes between the oceanic and global solutions 
to degree 10 we observe a higher energy in the oceanic solution, which was 
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expected, but also a similar decay in energy with increasing degree?. Such a 

behavior and the fact that many of the SST features above degree 10 have 

been filtered out by Levitus suggest that the spectral content of the SST 
could be approximated by the estimates of Table 1 reasonably well. 

The three sets of coefficients given in the Tables as well as the 

coefficients of the global SET3 solution to 4 max = 36, have been used to 

generate SST estimates on a l*xl° grid. Contour maps of these estimates 
based on a 5*x5* grid were then generated and are shown in Figures 6-9 of 
Appendix A, while a contour map of the Levitus set (SET3) also based on a 
5*x5* is shown in Figure 5. Examining Figures 6 and 7, that portray the SST 
from the ocean and global solutions up to degree 10 respectively, it is seen 
that the primary difference is that the former is able to show the broad 

features of the Gulf stream in a much better way than the latter. Examining 
Figures 8 and 9 that portray the SST features to degree 36 we can see that 
all features of the SST are almost identical and have an excellent agreement to 
the ones of Figure 5. In a more detailed visual examination the SST in the 
Gulf stream region from the SET3 data set and the four solutions has been 
plotted and' is shown in Figures 10 through 14 of Appendix. Now one can 
clearly see the inadequacy of the global solutions to reproduce the SST 
features faithfully, although the global solution to degree 36 provides a good 
qualitative approximation. The ocean only solution up to degree 36 is identical 
to the original data. 

In order to create Figures 8 and 13 the SST values on land were replaced 
by the corresponding values of the global solution to 4, nax = 36, since the land 
values of the ocean solution to 4 max = 36 are on the order of 20 meters and 
no contouring could be made. The same values have been used to plot the 
Levitus maps in Figures 5 and 10. Such a substitution may have created 
small distortions close to the coastlines which though are expected to be 
minor. In any case the values on land and close to the coastlines are 
meaningless. 


Conclusions. 


From the analysis made so far it is understood that the determination of 
spherical harmonics from the oceanic Levitus data set is very sensitive to a 
number of factors. It is seen that high correlations exist between the 
recovered coefficients that makes them completely unreliable, although they 
provide a good fit to the data. The best set of coefficients among solutions 
has been identified (SET3 ocean solution to & max = 10) based on certain 

criteria. Additional solutions which are as good (SET1 and SET2 ocean 
solutions) have been found but have not been chosen because they also 
reflect the data in the north polar region and therefore they do not provide 
an equally good fit to the ocean areas also attained by altimetry, as the SET3 
solution does. The experience obtained from this investigation indicates that 
although the chosen solution is the best among solutions that have been 
attempted it may very well turn to be suboptimum. As a matter of fact, it 
may be possible that some other slightly different configuration of the data 
and/or solution may very well give better results. Based on the above, one 
could even be inclined to use the coefficients computed from the global 
solutions since such solutions although, not giving very good fits, have very 
well established properties. 
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Appendix A 

Sea Surface Topography Maps 
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Figure 6. SST Estimates from the Oeeiari Analysis of Levitus Data (SET3) up to Harmonic Degree 10, Based 

on a 5"x5 p Grid. C.I. = 20 cm. 
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Figure 7. SST Estimates from the Global Analysis of Levitus Data (SETS) up to Harmonic Degree 10, Based 

on a 5“x5” Grid- C.I. = 20 cm. 
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Figure 12. SST Estimates from the Global Analysis of Levitus Data (SET3) up to Harmonic Degree 10 in the 

Gulf Stream Region, Baaed on a l*xl* Grid. C.I. = 10 cm. 
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Appendix B 

Spherical Harmonic Coefficients from Different Solutions of the Levitus SST 

Estimates. Units are in meters. 


27 


TABLE 


1 



OCEAN SOLUTION BASED 

ON 

LEV ITUS SST 

(SET3 ) 
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DEGREE OF EXPANSION =10 



NUMBER OF 1 

JNKNOWNS 
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-05 
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OF 
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OF 
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TO 
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N 

M 
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3 
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0 
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0.0000 
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0.0368 


4 

3 


0.0010 
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1 
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-0.0416 


6 

5 


0.0283 
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7 
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8 
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0.0102 



SN(CUM) 
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8 
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8 
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8 
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TABLE 2 


GLOBAL SOLUTION BASED ON LEVITUS SST (SET3 ) 


MAXIMUM DEGREE OF EXPANSION 
NUMBER OF UNKNOWNS 
RECIPROCAL CONDITION NUMBER 
DETERMINANT 


10 

121 

0 . 948E-01 
3.353*10.0**575 


NUMBER OF SST VALUES USED = 64800 
DEGREES OF FREEDOM - 64679 
RMS FIT TO LEVITUS DATA = 0.1 9M 


N M 

0 0 

1 0 

1 1 

2 0 

2 1 

2 2 

3 0 
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4 4 

5 0 
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5 3 

5 4 

5 5 

6 0 

6 1 
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SN 
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- 0.0076 
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0.0431 


0.0415 


0.3693 


0.3716 
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TABLE 3 


OCEAN SOLUTION BASED ON LEVITUS SST (SET3 ) 


MAXIMUM DEGREE OF EXPANSION 
NUMBER OF UNKNOWNS 
RECIPROCAL CONDITION NUMBER 
DETERMINANT 


36 

1369 

0 . 134E-08 
4.054*10.0**4635 


NUMBER OF SST VALUES USED * 301153 
DEGREES OF FREEDOM ■ 28784 

RMS FIT TO LEVITUS DATA * 0.01M 


N M 

0 0 

1 0 

1 1 

2 0 

2 1 

2 2 

3 0 

3 1 

3 2 

3 3 

4 0 

4 1 

4 2 

4 3 

4 4 

5 0 

5 1 

5 2 

5 3 

5 4 

5 5 

6 0 

6 1 

6 2 

6 3 

6 4 
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3 
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8 

4 
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8 

5 
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8 

6 
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8 

7 

- 0.3193 

8 

8 

- 0.2158 

9 
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- 0.2322 

9 

1 

- 0.3324 

9 

2 

0.3492 

9 

3 

0.1281 

9 

4 

0.1157 

9 

5 

- 0.0037 

9 

6 

0.2248 

9 

7 

- 0.1672 

9 

8 

- 0.0331 

9 

9 

- 0.1213 

10 

0 

0.1158 

10 

1 

- 0.2882 

10 

2 

- 0.1185 

10 

3 

- 0.1696 

10 

4 

0.1137 

10 

5 

- 0.0517 

10 

6 

0.2827 

10 

7 

- 0.0207 

10 

8 

0.0746 

10 

9 

0.1409 

10 

10 

- 0.0488 

11 

0 

0.1773 

11 

1 

- 0.0220 

11 

2 

- 0.2388 

11 

3 

- 0.2535 

11 

4 

- 0.0712 

11 

5 

- 0.0541 

11 

6 

0.2354 

11 

7 

0.1608 

11 

8 

- 0.0430 

11 

9 

0.2526 

11 

10 

- 0.0344 

11 

11 

0.0947 

12 

0 

- 0.2445 

12 

1 

0.3974 

12 

2 

- 0.0178 

12 

3 

- 0.1595 

12 

4 

- 0.0821 

12 

5 

- 0.0276 

12 

6 

0.0632 

12 

7 

0.4752 

12 

8 

0.0126 

12 

9 

0.0874 

12 

10 

- 0.1428 

12 

11 

- 0.0543 

12 

12 

0.1495 

13 

0 

- 0.5600 

13 

1 

- 0.1256 

13 

2 

0.3153 

13 

3 

0.2417 

13 

4 

0.1825 

13 

5 

0.0708 

13 

6 

- 0.2434 

13 

7 

0.3933 


0.0998 

- 0.5187 

- 0.4279 

0.0079 

0.1606 

- 0.0902 

0.0000 1.0952 

- 0.0399 
0.0838 
0.4949 
- 0.3580 
- 0.5554 
- 0.2888 
0.0939 
0.0494 
0.0640 

0.0000 1.0206 

- 0.0943 
- 0.3587 
0.5589 
0.0414 
- 0.4036 
- 0.3418 
- 0.1467 
- 0.0531 
0.0492 
0.1429 

0.0000 0.7927 

0.0910 
- 0.2856 
0.1834 
0.3351 
- 0.0478 
- 0.1418 
- 0.1409 
- 0.1739 
- 0.0254 
- 0.0248 
0.0255 

0.0000 0.8285 
0.0672 
- 0.1147 
- 0.1500 
0.1711 
0.2169 
0.0824 
0.0615 
- 0.1290 
- 0.0288 
- 0.0360 
- 0.0702 
- 0.0374 

0.0000 1.0991 

- 0.0220 
0.4347 
0.1112 
- 0.1929 
0.0587 
0.1462 
0.0389 


3.5653 


3.7085 


3.7923 


3.8818 


4.0344 
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13 

8 

0.3083 

0.0020 


13 

9 

0.0619 

0.0068 


13 

10 

- 0.1814 

- 0.0349 


13 

11 

- 0.1608 

- 0.1641 


13 

12 

0.0799 

0.0889 


13 

13 

- 0.0483 

- 0.0261 


14 

0 

- 0.1410 

0.0000 

1.4896 

14 

1 

- 0.7008 

- 0.1921 


14 

2 

0.2002 

0.1863 


14 

3 

0.2795 

0.7613 


14 

4 

0.4308 

0.0627 


14 

5 

0.3277 

- 0.2956 


14 

6 

- 0.2647 

- 0.0425 


14 

7 

- 0.1052 

- 0.2020 


14 

8 

0.2948 

- 0.0513 


14 

9 

0.2235 

0.1067 


14 

10 

- 0.1484 

0.0308 


14 

11 

- 0.0873 

- 0.3348 


14 

12 

0.2035 

0.0788 


14 

13 

- 0.0326 

0.2180 


14 

14 

- 0.2012 

- 0.0102 


15 

0 

0.7690 

0.0000 

1.8346 

15 

1 

- 0.6433 

- 0.1638 


15 

2 

- 0.2020 

- 0.2195 


15 

3 

0.2248 

0.5732 


15 

4 

0.4069 

0.6615 


15 

5 

0.3853 

- 0.2367 


15 

6 

0.0537 

- 0.3320 


15 

7 

- 0.2872 

- 0.3256 


15 

8 

- 0.0886 

- 0.3603 


15 

9 

0.1667 

0.1460 


15 

10 

- 0.0156 

0.2967 


15 

11 

- 0.1083 

- 0.2750 


15 

12 

0.4128 

- 0.2692 


15 

13 

0.0527 

0.2832 


15 

14 

- 0.2127 

0.0514 


15 

15 

- 0.1377 

0.0382 


16 

0 

0.9291 

0.0000 

1.8639 

16 

1 

0.3775 

0.0953 


16 

2 

- 0.3550 

- 0.6428 


16 

3 

- 0.1938 

- 0.3851 


16 

4 

0.2094 

0.5918 


16 

5 

0.0596 

0.3261 


16 

6 

0.1954 

- 0.3998 


16 

7 

0.0950 

- 0.2185 


16 

8 

- 0.1285 

- 0.4780 


16 

9 

- 0.1578 

- 0.0679 


16 

10 

0.1658 

0.4254 


16 

11 

- 0.2397 

0.1705 


16 

12 

0.3792 

- 0.3115 


16 

13 

0.2557 

0.0518 


16 

14 

- 0.2179 

0.0050 


16 

15 

- 0.1642 

- 0.1460 


16 

16 

- 0.0231 

0.0295 


17 

0 

0.1095 

0.0000 

1.8368 

17 

1 

0.8762 

0.1793 


17 

2 

0.0390 

0.0258 


17 

3 

- 0.2581 

- 1.0306 


17 

4 

0.1009 

- 0.4557 


17 

5 

- 0.2689 

0.4943 



4.3006 


4.6755 


5.0334 


5.3581 
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17 

6 

- 0.2565 

17 

7 

0.2932 

17 

8 

0.3655 

17 

9 

- 0.2050 

17 

10 

0.1640 

17 

11 

- 0.0972 

17 

12 

- 0.0107 

17 

13 

0.3555 

17 

14 

- 0.1302 

17 

15 

- 0.1085 

17 

16 

- 0.1207 

17 

17 

0.0193 

18 

0 
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18 

1 

0.3064 

18 

2 

0.3324 

18 

3 

- 0.2019 

18 

4 

- 0.0111 

18 

5 

- 0.0045 

18 

6 

- 0.7297 

18 

7 

- 0.1099 

18 

8 

0.5869 

18 

9 

0.1866 

18 

10 

0.0412 

18 

11 

0.2244 

18 

12 

- 0.2985 

18 

13 

0.0776 

18 

14 

0.0177 

18 

15 

0.0972 

18 

16 

- 0.0402 

18 

17 

- 0.1542 

18 

18 

0.0037 

19 

0 

- 0.5070 

19 

1 

- 0.4179 

19 

2 

- 0.0305 

19 

3 

0.1781 

19 

4 

- 0.1571 
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0.0530 

34 

3 
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- 0.0031 
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0.0224 
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- 0.0093 
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- 0.0002 

0.0008 
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0.0075 
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0.0199 
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- 0.0044 

- 0.0030 
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0.0010 
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- 0.0023 
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- 0.0050 
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- 0.0024 
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